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, Abstract. We study the fifth term in the asymptotic expansion of the heat operator 

• trace for a partial differential operator of Laplace type on a compact Riemannian 
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manifold with Dirichlet or Neumann boundary conditions 



il Statement of results 



D = -(Tr(V 2 ) + E). 



Let M be a smooth compact Riemannian manifold of dimension m with smooth 
boundary dM. Let V be a smooth vector bundle over M and let D be a partial 
J> . differential operator of Laplace type on C°°(V). There exists a unique connection 

0^ ' V on V and a unique endomorphism E of V so that 

(N ' 
O 

o 

' If 4> € C°°(y), let </> ;m be the covariant derivative of 4> with respect to the inward 

Q\ , unit normal. Let S be an endomorphism of V'jaM- We define the modified Neumann 

boundary operator B$ and the Dirichlet boundary operator B~ by 

Q_i' Bg4> := (4>-, m + S(p)\dM and B~4> := 4>\dM- 

Let B = Bg denote either boundary operator; we set S = for Dirichlet boundary 
conditions to have a uniform notation. Let D& be the operator defined by the 
appropriate boundary condition. If F is a smooth function on M, there is an 
asymptotic series as t J. of the form 

Tr L *(Fe- tD °) S S„> t ( ^ m)/2 a„(^ D,B) 

where the a n (F,D,B) are locally computable; see [6] for details. We computed a n 
for n < 4 in [3]; we have changed notation slightly from that paper. In this paper, 
we compute 0,5 if the boundary is totally geodesic. If the boundary is not totally 
geodesic, the second fundamental form L enters and the combinatorial complexity 
becomes formidable. In this setting, we also obtain a result by restricting attention 
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to regions in R m . In particular, we determine the coefficient of |VL| 2 in ^(Aq, Bq) 
for the Laplacian A on functions with pure Neumann or Dirichlet boundary con- 
ditions. In the Neumann case, this coefficient plays an important role in the work 
of Gursky [7] on compactness of isospectral sets in planar domains. The determi- 
nation of this coefficient was, in large part, the original motivation for the present 
paper. 

We adopt the following notational conventions. Let indices i, j, ... range from 1 
through fii and index a local orthonormal frame {ei, e m } for the tangent bundle 
of M. On the boundary, let indices a, b, ... range from 1 through m — 1 and index 
a local orthonormal frame for the tangent bundle of the boundary; let e m be the 
inward unit normal. We adopt the Einstein convention and sum over repeated 
indices. Let Rijki be the components of the curvature tensor of the Levi-Civita 
connection on M. With our sign convention, #1212 is negative on the standard 
sphere in Euclidean space. The Ricci tensor p and the scalar curvature r are given 

by 

Pij '■= Rikkj and T = pa = Rikki- 

Let p 2 := PijPij and R 2 '■— RijkiRijki be the norm of the Ricci and full curvature 
tensors. Let Clij be the endomorphism valued components of the curvature of the 
connection on V, and let L a b be the components of the second fundamental form; 

L a b '■— (V ea e b , e m ). 

Let ';' denote multiple covariant differentiation with respect to the Levi-Civita 
connection of M and let ':' denote multiple tangential covariant differentiation on 
the boundary with respect to the Levi-Civita connection of dM; the difference 
between these two is measured by the second fundamental form. When sections 
of bundles built from V are involved, ';' will mean V M <g> 1 + 1 <g> V y and ':' will 
mean V aM <£> 1 + 1 <£> V y ' 3M . Thus E. a — E :a since there are no tangential indices 
in E to be differentiated, while E. ao and E :ao do not agree since the index a is also 
being differentiated. Since L and S are only defined on the boundary, these tensors 
can only be differentiated tangentially. Let dx and dy be the Riemannian volume 
elements on M and on dM. If /1 G C°°(M) and if f 2 e C°°(<9M), define 

h[M] := J M h(x)dx and f 2 [OM] = J dM .h(y)dy. 

Theorem 1.1. Suppose the boundary of M is totally geodesic. Then 

a 5 (F,D,B^) = ±5760- 1 (47r) (m_1)/2 Tr{F(360£; ;mm + U40E [m S + 720E 2 + 
2A0E :aa + 240tE + 120n ab n ab + 48r ;l4 + 20r 2 - 8p 2 + 8R 2 - 120p mm E - 2Qp mm r + 
480rS' 2 + (90+, -360-)n am n am + 12r ;mro + 24/? 

mm:aa H~ ^^Pmm;mm 

120p mm S + 960S :aa S + 600S :a S :a + 16RammbPab — 17 PmmPmm ~ WRammbRammb + 

2880.BS' 2 + 1440S' 4 ) + J F ;m ((195/2+, 60")T ;m + 240t5 - 90p rnm S + 270S :aa + 
(630+,450-)£; ;m + 1440£S' + 7205 3 ) + F ;mm (60r - 90p mm + 360E + 360S 2 ) + 
180SF ;mmm + 4:5F ;mmmm }[dM}. 

As mentioned above, if dM is not totally geodesic, the number of invariants be- 
comes unmanageable; the general formula involves 63 additional terms; see Lemma 
5.1 for details. However, if we restrict to the Laplacian Ao on functions with Neu- 
mann or Dirichlet boundary conditions, assume M is a domain in R m , and set 
F = 1, a great simplification occurs. Recall that L oc . a — L ac - = R a bcm, see for 
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example [3, Lemma A.l (b)]. Consequently, if the curvature tensor R vanishes, 
then VL is a symmetric 3-tensor. Thus there exist universal constants such that 
a 5 (A ,B±) = ±(4Tr)( m - i y 2 5760- 1 J{e^L ab:c L ab:c + efL aa L bb L cc L dd 

+ e 3 LabLatLccLdd + e ^ L ab L ab L c dL c d + L ab L bc L ca L dd 
+efL ab L bc L cd L d a}[dM}. 

Theorem 1.2. e+ = -19 • 45/16 and e{ = 45/16. 

Remark:. The constant e\ is of particular interest as it controls certain compact- 
ness estimates for Neumann boundary conditions and plays an important role in 
theorems of Gursky [7]; it was not previously known and its determination was 
one of the primary goals of this study. The value of was previously known; see 
[11,12]. Levitin [10] has used the algorithm of Kennedy [8, 9] to compute a$ for the 
ball in R m for any to; this result has been used by van den Berg [1] to determine 
e 2 , e*, (a + es)*, and ef . See also related work by Bordag et al [2]. 

Here is a brief outline to the paper. In §2, we recall some functorial properties of 
the invariants a„(-), and recall the calculation of a n for n < 4. In §3, we compute a§ 
for the special case m = 1; this is an important step in the general case. In §4, we 
prove Theorem 1.1; we shall omit the calculation of the coefficient of O am f2 am as 
it is fairly lengthy; details arc available from the authors. In §5, we prove Theorem 
1.2. §6 is an appendix in which we give some variational formulas used elsewhere 
in the paper. 

§2 Functorial properties 
We summarize below properties established in [3] : 
Lemma 2.1. 

1) Let N V (F) — F ;m ... m be the v th normal covariant derivative. There exist invari- 
ant local formulae a n (x,D) and a n ,u(y,D) so that: 

a n (F,D,B) = {Fa n (x,D)}[M] + {Z < v < n _ 1 N v {F)a n , v {y,D,B)}[dM]. 

The interior invariant a n {x,D) is homogeneous of order n in the jets of the 
symbol of D and vanishes for n odd; it is independent of the boundary condition 
chosen. The boundary invariants a nj „ are homogeneous of order n — v — 1. 

2) Let M = Mi x M 2 and D = D 1 <g> 1 '+ 1 <g> D 2 . Let dM x = 0. Then 

a n (x,D) = H p+q=n a p {xi, D 1 )a g (x 2 , D 2 ) 

a n ,v(y,D) = Z p+ q =n ap tL ,(x 1 ,D 1 )a q (y2,D 2 ,B). 

3) If we expand a n (x,D) or a n ^(y,D,B) with respect to a Weyl basis, the coeffi- 
cients depend on the dimension to only through a normalizing constant. They 
are independent of the dimension of V. 

4) IfD(e) = e- 2 * F D, then £\ e=0 a n {l, D,B) = (to- n)a n {F,D,B). 

5) IfD(e) = D-eF-I v , then £\ e=0 a n (l, D,B) = a n . 2 (F, D,B). 

6) IfD(e) = D-e-I v , then £\ e=0 a n (F,D,B) = a n . 2 (F,D,B). 

7) Ifm = n + 2, then £ \ e=0 a n {e- 2e f F, e~ 2,if D,B) = 0. 



4 



T. BRANSON, P. GILKEY, AND D. VASSILEVICH 



8) Let m = 1 and let b G C°° [0, 1] be real. Let A = d x -b, A* = -d x -b, D x = A* A, 
D 2 = AA*, and S = b. Then 

(n - l)(a n (F, D\,B~) - a n (F, D 2 ,B+)) = a n _ 2 (d 2 x F + 2bd x F, D U B~). 

Remark. Assertion (3) is a bit formal; it is best illustrated by reference to Theo- 
rems 2.4 and 2.5 below which express a n for n < 4 in terms of a Weyl basis. 

There are two additional properties which we shall need which were not dis- 
cussed in [3]. The first relates certain coefficients in a Weyl basis for Neumann 
boundary conditions to the corresponding coefficients for Dirichlct boundary con- 
ditions. Suppose that M = N x S 1 and that V is the trivial line bundle. Let 
T{y,0) = (y,—0) define an involution of M where 9 £ M./2-k'Z — S 1 is the usual 
periodic parameter. Let M = N x [0, n] C M and let F be a smooth function on 
M which is preserved by T. Let e m be the inward unit normal on M; e m = d$ at 
8 = and e m — — de at 9 — n. We assume D is preserved by T; this means that if 
D = -{g^didj + A k d k + B], then 

g ab (y, 9) = g ab (y, -9), g mm (y, 9) = g mm (y, -9), 

g am (y,9) = -g am (y,-9), A m (y,9) = -A m (y,-9), 

A a (y, 9) = A a (y, -0), B(y, 9) = B(y, -9). 

Let denote pure Neumann and Dirichlet boundary conditions so that S = 0. 

Lemma 2.2. We adopt the notation established above. Let n be arbitrary and let 
v be even. Then 

a n Ay> D >®o) + an,v(y,D,B~) = 0. 

Proof. Let F be an even function. Since D is invariant under the involution T, we 
may decompose the eigenfunctions of D on M into the even and odd eigenfunctions. 
The even eigenfunctions satisfy pure Neumann boundary conditions on M; the odd 
eigenfunctions satisfy Dirichlct boundary conditions on M. Thus 

Tr L 2 M (Fe- tD -) +Tr L 2 M (Fe- tD +) = Tr L2[£{) (Fe- tD ). 

We equate coefficients of t in the asymptotic expansions of both sides of this equa- 
tion. The interior integrals cancel since on the left hand side we are integrating 
over [0,7r] twice and on the right hand side we are integrating over [— n, n]. There 
are no boundary integrals on the right hand side so the boundary integrals on the 
left hand side must cancel. This shows 

Z„{N»(F)(a n Ay,D,B+) +a n , v {y,D,B-))}[dM] = 0. 

Since we may specify the N V (F) arbitrarily for v even, the desired vanishing theo- 
rem now follows. □ 

In the next Lemma, we shall restrict to the case n = 5 in the interests of clarity 
and note that there is a more general principal applicable. 
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Lemma 2.3. When as is expanded in a Weyl basis, the coefficients of the follow- 
ing terms are zero: F- m Tr(f2 am;a ), F Tr(ri am;ma ), FTr(SQ am;a ), FTr(S :a n am ), 
FL aa Tr(Q bm;b ). 

Proof. We adapt an argument given in [4] . Let F be real. If D and S are real, then 
Tr L 2(Fe~* DB ) is real so the coefficients of the terms listed above must be real. On 
the other hand, if V is unitary with respect to some fiber metric on V and if E and 
S are self-adjoint, then is sclf-adjoint so again Tr L 2(Fe~ tDB ) is real. We take 
V to be a line bundle; f2 is pure imaginary in this context. Thus the coefficients 
must be pure imaginary as well and hence must vanish. □ 

We refer to [3] for the computation of a n for n < 4. 

Theorem 2.4 (Dirichlet). 

1) a (F,D,B-) = (4 7 r)- m / 2 Tr( J F 1 )[M]. 

2) ai (F,D,B-) = -4- 1 (4 7 r)-( m - 1 )/ 2 Tr(F)[9M]. 

3) a 2 {F,D,B-) = (47r)~" l//2 6~ 1 {Tr(6F_E + Ft) [M] + Tr(2FL aa — 3F- m )[dM]} 

4) a 3 {F,D,B-) = -384- 1 (47r)-( m - 1 )/ 2 {Tr{F(96S + 16r - 8p mm + 7L aa L bb - 
\{)L ab L ab ) - 30F. m L aa + 24F, }[dM]} 

5) a 4 {F,D,B-) = (47r)- m / 2 360- 1 {Tr{f 1 (60 J B :fefc + 60T£;+180 J B 2 + 30f7 2 + 12T :fefc + 
5T 2 -2p 2 +2i? 2 )}[Af]+Tr{F(-120F ;m -18r ;m +120SL aa +20rL Qa +4i? amam L fc;) - 
l2R ambm L ab + 4:R abcb L ac + 24:L aa:bb + 0L ab:ab + AQ/2\L aa L bb L cc + 0f2 am;a — 
88/7L ab L ab L cc +320/21L ab L bc L ac )+F. m {-l80E-30T+0R amam -180/7L aa L bb + 
60/7L ab L ab ) + 24F, mm L aa )}[dM}} 

Theorem 2.5 (Neumann). 

1) a (F,D,B+) = (47r)-™/ 2 Tr(F)[M]. 

2) ai {F,D,B+) = 4- 1 (47r)-( m - 1 )/ 2 Tr(F)[<9M]. 

3) a 2 {F,D,B+) = (4 7 r)- m / 2 6- 1 {Tr(6FS + FT)[M]+Tr(2FL aa + 3F :m + 12F5) 

[dM]}. 

4) a 3 (F,D,B+) = 384- 1 (4^)-(" l - 1 )/ 2 {Tr{F(96£ + 16r - 8p mm + l3L aa L bb + 
2L ab L ab + 96SL aa + 192S 2 ) + F ;m (6L aa + 96S) + 2AF. mm }[M]} 

5) a 4 {F,D,B+) = (4 7 r)- m / 2 360- 1 {Tr{F(60£; :fefe + 60rS+180 J B 2 + 30f7 2 + 12T ;fefe + 
5r 2 - 2p 2 + 2i? 2 )}[M] + Tr{F(240£ ;m + 42r ;ro + 24L aa:bh + 0L ab:ab + 120EL aa + 
20t L aa +AR amam L bb — 12R ambm L ab + 4R abcb L ac + 40/3L aa L bb L cc + 8L ab L ab L cc + 
32/3L ab L bc L ac +0Q am . a + 720SE+120ST+0SR amam + lUSL aa L bb +48SL ab L ab 
+A80S 2 L aa + A80S 3 + 120S :aa ) + F :m (l80E + 30t + 0R amam + \2L aa L bb + 
12L afc L Qb + 72SL aa + 2AQS 2 ) + F. mm (24L aa + 1205) + 30F. llm }[dM] } 

§3 Dimension One 

This section is devoted to the calculation of a$ in the one dimensional case; this 
forms the foundation for our later computations. We restrict to scalar operators; 
recall S = if B = B~. 

Theorem 3.1. a 5 (F,D,B^) = ±384- 1 {(24£; :mm + 96E, m S + 48E 2 + 192ES 2 + 
96S 4 )F + ((42+, 30-)£7 ;m + 96ES + A8S 3 )F :m + (24E + 2AS 2 )F :mm + 12SF. mmm + 
3F. }[dM}. 
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Remark. If we set F = 1, this agrees with Theorem 3.4 of [5]. 

Proof, we use Lemma 2.1 (1) to see there exist universal constants so that 
a 5 (F, D, Bf) = ±38A- 1 {(bfE. mm + b+E, m S + bfE 2 + b+ES 2 
+b+S+)F + {bfE, m + b+ES + 6+S 3 )F ;m + (bfE + b+ S 2 )F., mm 
-\-b+iSF; mmTn ^2^;™™}!^]' 
Let M be the interval [0, n], and set S = 0. Double the interval to obtain a copy 
of the circle, and assume E and F are even on the double. We apply Lemma 2.2 
to see 

a 5 (F,D,B+)+a 5 (F,D,B-) = 0. 

This shows that b + = 6^, b+ = b^ , b+ = bg , and b+ 2 = b^ 2 . No conclusion can be 
drawn about the terms involving S since we have set S = 0. Also, no conclusion can 
be drawn about the term involving E- m since this will vanish for an even function. 
This shows that the general formula for as in the one-dimensional case is 

a 5 (F, D, B+) = ±384- 1 {(b 1 E. mm + b+E, m S + b 3 E 2 + b+ES 2 

+b+S 4 )F + (b+E. m + b+ES + b+S 3 )F. m + (b 9 E + b+ a S 2 )F [mm 

+b+ 1 SF- immm + bi2F. mmmm }[dM}. 

By Theorems 2.4 and 2.5, 

a 3 (F, D, B+) = ±384- 1 {(96£ + 192S 2 )F + 96SF :m + 24F ;mro }[3M]. 

Let D{e) = D-el. Then E{e) = E + el. By Lemma 2.1 (6) 

f e \ eM (F,D,B+)=a 3 (F,D,Bf). 

This shows b 3 = 48, b+ = 192, b+ = 96, and b 9 = 24. Next let D(e) = D — eF. Then 
E(e) = E + eF. By Lemma 2.1 (5) 

£\ e=0 a 5 (l,D,B±) = a 3 (F,D,B+). 

It now follows that b\ = 24 and b+ = 96. This shows that 

a 5 {F, D, Bf) = ±384- 1 {(24S ;mm + 96E- m S + A8E 2 + 192ES 2 + b+S A )F 

+ (b±-E. m + 96ES + b+S 3 )F. m + (24£ + b+ S 2 )F. mm 

+b 11 SF- mmm + b^F. 

mmmm 

}[dM]. 

We now adopt the notation of Lemma 2.1 (8). Let M = [0, 1] and let b e C°°[0, 1] 
be real. Let A = d x — b, D\ = A* A, and D 2 = AA*. Then E 1 = -b x - b 2 and 
F<2. =b x — b 2 . We take Dirichlet boundary conditions for D\ and modified Neumann 
boundary conditions given by A* for D 2 ; this means S = b. We assume F vanishes 
to infinite order at the endpoint x = 1. We expand a 3 with Dirichlet boundary 
conditions, 

a 3 (F xx + 2bF x ,D 1 ,B-) 

= -384- 1 {96(-b x - b 2 ){F xx + 2bF x ) + 2A(F XX + 2bF x ) xx }(Q) 

= -384- 1 {(-966 :c f; ;c - 9&b 2 F xx - l92b x bF x - 1926^ 

+ 96b x F xx + 48bF xxx }(0) 

= -38A- 1 {2AF XXXX + A8bF xxx - 96b 2 F xx 

+ (48b xx - 1926,6 - 1926 3 )F X }(0), 
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a$ with Neumann boundary conditions, 

-a 5 (F, D 2 ,B+) = -384" 1 {(24(6 X - b 2 ) xx + 96(6, - b 2 ) x b 
+48(6, - 6 2 ) 2 + 192(6, - 6 2 )6 2 + 6+6 4 )F 
+(b+(b x - b 2 ) x + 96(6, - 6 2 )6 + b+b 3 )F x 
+ (24(6, - 6 2 ) + b+ b 2 )F xx + 6+6F„, + 
and as with Dirichlet boundary conditions, 

a 5 (F, D,B~) = -384- 1 {(24(-6, - 6 2 )„ + 48(-6, - 6 2 ) 2 )F 
+6-(-6, - 6 2 ),,F, + 24(-6, - b 2 )F xx + 6 12 F,,,,}(0). 
By Lemma 2.1 (8), 

a 3 (F xx +2bF x ,D u B-) = 4{a 5 (F,D u B-) - a 5 (F, D 2l B+)}. 
This leads to the equations 



Term 


Equation 


Term 


Equation 


b 4 F 


= 48 - 192 + bt + 48 


6 3 F, 


-192 = 4(-96 + 6+) 


b x bF x 


-192 = 4(-26+ + 96 - 26 fi ~) 


6 F 

U XX ± X 


48 = 4(6+ - 6 fi ") 


b F xx 


-96 = 4(-24 + 6+ - 24) 


hF 

u± XXX 


48 = 46+ 


± xxxx 


24 = 4(6 12 + 6 12 ) 







We solve this system of equations and complete the proof by computing: 

6+ = 96, bt = 42, 6g = 30, 6+ = 48, 
bf = 24, 6+ = 12, 6 12 =3. □ 

§4 Totally geodesic boundary 

This section is devoted to the proof of Theorem 1.1. We assume for the remainder 
of this section that the boundary of M is totally geodesic. We begin the proof of 
Theorem 1.1 by expressing as in terms of universal expressions with undetermined 
coefficients. Note that 5760 = 15 • 384 = 16 • 360. 

Lemma 4.1. Suppose the boundary of M is totally geodesic. There exist universal 
constants so that a 5 (F, D, £+) = ±(4 7 r)( m - 1 )/ 2 5760- 1 Tr{F(360 J & mm +1440£ ;m S+ 
720E 2 + 2880£S 2 + 14405 4 + 240£ :aa + 240t£ + I20n ab n ab + 48r ;4l + 20r 2 - 
8 P 2 + 8i? 2 - 120p mro E - 20/wr + 480t5 2 + cfn am n am + 

C-&Pmm\mm H~ C5 Ty m S ~\~Cq Pmm$ H~ C7 S- aa S -\- Cg S- a S- a H~ C9 Rammb Pab H~ ClO PmmPmm H~ 

CllRammbRammb + F. m ( (630+ , 450~ )E, m + 1U0ES + 7205 3 + cf 2 T, m + 2A0tS + 
ctsPmmS + c+ 4 S :oa ) + F ;mm (360£ + 360S 2 + 60r + c 15 p mm ) + 18QSF ;mmm + 

Proof. Since the boundary is totally geodesic, we may interchange ';' and ':' as 
convenient. We write down a basis for the space of invariants using H. Weyl's 
theorem. The only tricky part of this analysis is the study of terms of order 4 in 
the jets of the metric; these terms are studied at the end of the proof. We omit the 
curvature terms R ma bcRmabc, and p mc Pmc since these vanish on a totally geodesic 
boundary. Since r ;m = 2p mm - mi we also omit the terms p m m-mS and F- m p mm]m 
from our list. Finally, we use Lemma 2.3 to omit the following invariants from our 
list: 

5"), Ti(Fn am S- a ), Tr(F ;m fi am;a ). 
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We can use previous results to determine many of the coefficients in a Weyl 
spanning set: 

1. We use Theorem 3.f to determine the coefficients of the following expressions 
which do not involve f2, R, or the tangential derivatives of S and E : 
FE ;mm , FE- m S, FE 2 , FES 2 , FS 4 , F- m E [m , F- m ES, 

J3 C 1 TP TP c2 



E-m,S j F- m7n E , F-mjnS , F- mmm S 1 F-^ 



\mmmm • 



2. We use Lemma 2.1 (6) to see j^\ t=0 a${l,D — eI,B) = 03(1, D,B); this enables 
us to determine the coefficient of tE and p mm E in as from the expression for 03 
given in Theorems 2.4 and 2.5. 

3. We apply Lemma 2.1 (2) to a product M = M x x M 2 where dM x = to see 

a 5 ,u(x 1 ,y 2 ,D 1 + D 2 ,Bg ) = T, p+q=5 a p (x 1 ,D 1 )a q ^(y 2 ,D 2 ,Bg). 

(i) We take (mi 7 m 2 ) — (4, 1) and study the terms in 05,0 which are independent 
of M 2 . These terms arise when (p, q) = (4, 1) so 

a 6 ,o(-,B| ) - ±(4 7 r)(™- 1 )/ 2 {60^ + 30(f^) 2 + 12r^ 

+ 5(t Mi ) 2 - 2(p Ml ) 2 + 2(i? Ml ) 2 + ...}/(4 • 384). 

This determines the coefficients of the invariants E- aai f2 a {,f2 a fc, t-u, t 2 , p 2 , 
and R 2 in 015,0; still to be determined, of course, is the coefficient of T ;mm . 
Similarly, we study the cross terms to evaluate the coefficient of tS 2 in a 5i o, 
to evaluate the coefficient of tS in 05,1, and to evaluate the coefficient of r in 
05^2; the cross terms arise when (p, g) = (2, 3). 

(ii) We take (mi,m 2 ) = (2,2) and compute the cross terms involving the curva- 
ture tensor in 05,0- These arise when (p, q) = (2, 3) so 

a 5l0 (-,B|) - ±(4 7 r)( m - 1 )/ 2 T Ml (16r M2 - 8p^)/(6 • 384) + .... 

Since r 2 = 2t Mi t M2 + ... and rp„ lm = r Ml p^ 2 m + the coefficient of t 2 
in a 5 is 16 • 5760/(2 • 6 • 384) = 20 and the coefficient of rp mm in 05.0 is 
-8 • 5760/(6 • 384) = -20. 
We use Lemma 2.2 to see that c+ = with our normalizing sign convention for 
many of the remaining coefficients; a notable exception is the coefficient of flam^am 
since this term vanishes under the hypothesis of Lemma 2.2. 

We conclude the proof of Lemma 3.1 by studying the terms which are linear in 
the second covariant derivative V 2 R of the curvature tensor. We must set some of 
the indices equal to m and contract the remaining indices in pairs in an expression 
Ru i 2 i3u;*5i6- If none °f t ne indices are equal to to, we must contract 3 pairs of 
indices. The usual argument then shows this result is a multiple of t-u. Thus we can 
restrict our attention to the case where either 2 or 4 of the indices i M are equal to 
to. If 4 of the indices are equal to m, we get a multiple of p TOTO;mm . In the remaining 
case, 2 of the indices are to and modulo lower order terms we must consider the 
invariants 

Rabba;mmi Ramma-.bbi Rabam;bmi and Rambm-ab- 

By the second Bianchi identity, = R a mbm;ab+Ramma;bb+Ramab;mb; this eliminates 
Rambm-ab from our list. Modulo lower order terms, we may replace R am ab;mb by 
Ramab-bm- By the second Bianchi identity, 

Rabam\bm ~t~ Rabmb;am ~t~ Rabba:mm — 2R a f Jam . f l )Trl -\- Rabba;mmi 
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this eliminate R a bam;bm from our list as well. We then note 

Pmm-.aai Pmm;mm} • ^ 

We complete the proof of Theorem 1.1 by evaluating the constants which appear 
in Lemma 4.1. We will not compute cf as this requires somewhat different tech- 
niques; details are available from the authors upon request. This term vanishes for 
the scalar Laplacian in any event. 

Let D(e) = e~ 2ef D, S(e) = e~ €f S, F(e) = e~ 2ei 1 F. We assume f- m \ dM = to 
preserve the condition that the boundary is totally geodesic. If to = 7, we use 
Lemma 2.1 (7) to see 

£\ eM (F,D) = 0. 

We use Lemma 4.1 and results from §6 to compute ^| £= oas(F, D). We integrate 
by parts and combine terms to express -^\ e= ga^{F, D) in terms of independent 
integrands and derive the system of equations 



Term 


Equation Term Equation 


F f;mmmm 


39 = 2c 2 + c 4 


F j *-,mmS 


120 = 4 


F f-.rnmaa 


303 = 6c 3 + 12c 2 + c 4 


F j - t mmT 


40 = 2c 2 + c g 


F f:aaPmm 


60 = -2c 3 - Cg - 2cio - C15 


F.f;aT;a 


12 = c 2 


F j \aPmm\a 


324 = -ci 5 + 12c 2 + 6c 4 


F\m J \aS -.a 


270 = c+ 4 


F f:aRammb-,b 


63 = 4c2 + c 4 


Ff-abPab 


16 = Cg 


F-mf-.mmm 


(195 + ,120") = 2c^ 2 


F;rnmf\aa 


90 = -ci5 


Ff-.aaS 2 


1080 = c+ + c+ 




-90 = c+ 


Ff;aS;aS 


360 = 4 - 4 


F f-mmmS 


270 = ct 


F f:abRammb 


378 = 5c 9 + 2cu + 24c 2 + 2c 4 


F f-.mmPmm 


346 = — ci 5 — 6ci — en + 12c 2 



We solve this system of equations and complete the proof by computing: 

c 2 = 12, c 3 = 24, c 4 = 15, c+ = 270, c+ = 120, 

4= 960, ct = 600, c 9 = 16, cio = -17, en = -10, 

c i2 = !95/2, Cl2 = 60, 4s = -90, c+ 4 = 270, c 15 = -90. □ 



§5 Domains in flat space 

In this section, we complete the proof of Theorem 1.2. We drop the assumption 
that the boundary of M is totally geodesic; thus ":" and ";" differ. Lemma 4.1 
generalizes to this context to become 

Lemma 5.1. Suppose the boundary of M is smooth but not necessarily totally 
geodesic. 

1) There exist universal constants so that 

a 5 (F, D, B|) = ±5760- 1 (4 7 r)( m - 1 )/ 2 Tr{F{360£; ;mm + lU0E. m S + 720E 2 + 
240E :aa + 240tE + 120ft ab ft ab + 48t :4J + 20t 2 - 8p 2 + 8R 2 - l20p mm E - 20p mm r + 
480tS 2 + (90+, -360-)fi aTO f) am + 12r ;mro + 24p 

mm:aa H~ ^^Pmm;mm 



span{i? a6ba; 
=span{-r ;ii , t. 



mm-, 

rami 
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120/9 

run i + 

2880ES 2 + 1440S 4 } + F :m {(195/2+, 60-)r :m + 240x5 - 90p mm S + 270S :aa + 
(630+,450-)S ;m + 1440ES + 720S 3 } + F ;mm {60r - 90p mm + 360E + 360S 2 } + 
180iS^ mmm +45i^ mmmm +i^{c2-L L aa E :rn -\-d,2 L aa T :rn -\-d^ L a bR a7nm b :rn -\~d^ LaaS-.bb~\~ 

$5 L a bS :a b + L aa: bS : b + d^ L a b: a S;b + d% L aa: bbS + Lab:abS + d^Laa-bLcc-b + 
di±Lab:a-L C c:b ~l~ d^-^abia-^bcic d^L a b:cL ab:c H~ d^L abic-L 'ac:b H~ d^L aa: bbL cc -\- 

diQL a b :a bL cc + d^yL a b :ac L/bc H~ d^L aa: bc L bc + digLbc-.aa^bc 

+ U40+L aa SE+ 

d2QL aa S p mm + 24Q + L aa ST + d2 1 L a t,Spab + d22L a bSR m a,bm + (195 + , 105 - ) L aa LbbE + 

(30+, 150-)i afc i afc S + (195+/6, 105-/6)L aa L 6b r + (5+, 2h~)L ah L ah T 

-\-d2^L aa LbbPmm ~t~ d2^L a bL a bPmm H~ d2^L aa LbcPbc H~ d2§L aa LbcRmbcm d2 f jL a bL a cPbc 
-\~d2^L a bL a cRmbcm H~ ^29^ afr-^ cdRacbd H~ d^QL aa S^ + d^L aa LbbS + d^-^ab-^ab^ ■ + 

d^L aa LbbL cc S + d^L aa Lb c LbcS + d^L a bLb c L ac S + d^L aa LbbL CC L <jd+ 
d^L aa LbbL c dL cc i + d^LabLabLcdLcd + d 3 QL aa Lb c L cc iLdb + d 40 L a bLb c L c ,iL ( ], a }+ 
^; m {(90+ -450^)L Qa i; + 4iaa/*mm + (15+, -75-)L aa r + d+ 2 L aa S 2 + dt 3 L aa:bb + 
d^Lab.ab + d 45 L a bp a b + d 4 gL a bR ma bm + d 47 L aa LbbS + d^L a bL a bS + d 4 gL aa Lb b L cc + 
d^LaaLbcLbc + d 51 L a bLb c L ac \ + F imm {d^ 2 L aa S + d 53 L aa Lbb + d 54 L a bL a b}+ 
F ;mmm df 5 L aa }[dM}. 

2) d+3 = 151.875, d+ 4 = -11.25, and = -30. 

3) dj 3 = 39.375, dj 4 = -11.25, and d^ = -105. 

Proof. We write down a basis for the space of invariants using H. Weyl's theorem. 
We use Theorem 1.1 to evaluate the coefficients of the expressions not involving 
the second fundamental form i; we also Lemma 2.1 (2) to evaluate other coeffi- 
cients using suitable product formulas. The following invariants can be expressed 
in terms of invariants already appearing in the list above and arc therefore omitted: 

E L aa -bRcmcbj EL a b: a Rcmcbi E L a bRacbc;nii E p 7nrn - m L aaj F/) mm ; m 5, F-^p^^-^. We 

can also replace ';' by ':' in some expressions at the cost of introducing additional 
lower order terms involving the second fundamental form. Thus, for example, by 
Lemma 6.1, we may express 

Pmm;m — Rabba\ra}j2 — (^~;m 2R a bbm;a) /2 

~( T ;m — 2p am : a + ■■■)/2 = ("T ;TO — 2Lbb:aa + 2L a b :a b + 

This completes the proof of 1). 

To prove 2) and 3), we suppress the coefficients of F and F :mm and we suppress 
terms of length greater than 1 in the coefficients of F ;m . This permits us to express 

a 5 (F,D,B$) = ±5760- 1 (47r)( m - 1 )/ 2 Tr{F(*) + 

+F ;m {(195/2+, 60-) • r, m + 270S :aa + (630+, 450-)£ ;m 

-\-d\^F- m L aa -bb + d^ 4 F- m L a b:ab + •••} + f;mm(*) 

+F ;mmm (180S + dj^Laa) + 45F ;mmmm } [0M] . 
We have to be a bit careful since we have not used linearly independent mono- 
mials in the variational formulas in §6. Also, and equally importantly, we must 
integrate by parts in computing certain integral invariants. We apply Lemma 2.1 
(7) to see that if m — 7, then 

0=£U= a 5 (e- 2 *fF,e- 2 *fD,B± {e) ). 

The information in the tables below uses information from §3 and from the appendix 
in §6; the sum of each column multiplied by the appropriate entry and coefficient is 
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zero; this yields the equations necessary to compute df 3 , df 4 , and df 5 and complete 
the proof. □ 



Neumann Boundary Conditions m = 7 



Term 


Coeff 






f-.a^-./m^ab-.b 


f \mmmF^m 


f;aam-F-m 


F E 


630 











2.5 


2.5 




97.5 











-12 


-12 


F^\m S-.aa 


270 





2.5 


-5 





2.5 


F 

L ;mmmm 


45 


-14 


-4 


-8 


-9 


-4 




180 


2.5 
















4> 





-4 


12 





-6 


^-.m^ab-.ab 


d+ 





-2 


10 





-1 


F L 


4, 


-6 


-3 












Dirichlct Boundary Conditions m = 7 



Term 


Coeff 




f:aF- rn Lbb:a 


f;aF- m L a b:b 


f \irtmmF- m ra 


f;aamF- m 




450 











2.5 


2.5 




60 











-12 


-12 


F\rnmmm 


45 


-14 


-4 


-8 


-9 


-4 


F- t mLaa:bb 







-4 


12 





-6 




^44 





-2 


10 





-1 


F;mmmL a a 




-6 


-3 












We can now proceed with our study. We set m = 2, set F = 1, and suppress all 
terms of length greater than 2 to see 

a 5 (l,D,Bf) = ±5760- 1 (4^) (m - 1)/2 Tr{(135/2)r ;mm 
+270r ;m 5' - 360$ i Si + dfL n E im + efL 11:1 L U :i 
+e^L 11:11 S + efL n r + ...}. 



Lemma 5.2. e£ = -53.4375, = 270, e+ = 22.5, e{ = 2.8125, = -90. 

Proof. We consider the variation D(e) — e~ 2ei D with B s (e) defined appropriately 
and with E(0) = r(0) = 0. By Lemma 2.1 (4), 

-3a 5 (f,D,B±) = || e=0 a 5 (l,e- 2 ^£»,B| (e) ). 

We use this relation to evaluate d 2 , e 1 , and ey from the tables given below; the 
variational formulas used to create these tables are contained in §6. We also apply 
the results of §3 to evaluate certain of the coefficients. □ 

Neumann Boundary Conditions m = 2 
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Term 


Coeff 


j ;m J - / ll: 11 


,/ ;mmm J -'ll 


J ;m ' ;m 


J ;m»~':ll 


J ;mmmm 




-3a*,(f ■) 




-421.875 


90 


-292.5 


-810 


-135 


-540 




67.5 


-4 


2 


-4 





-2 







270 











-2 





-2 


^11:1-^11:1 


e i 


2 

















Ln-.uS 


e 7 











-1 








LllT- m 


e 8 


-2 


-2 


-1 












Dirichlct Boundary Conditions m = 2 



Term 


Coeff 


f;mLll:ll 


fymrnmL 11 


f;mT;m 


f;mmmm 


-3a 5 (/,-) 




-84.375 


315 


-180 


-135 


7"; mm 


67.5 


-4 


2 


-4 


-2 


^11:1^11:1 


er 


2 















-2 


-2 


-1 






§6 Appendix (Variational formulas) 

In the appendix, we present some identities used in the previous sections. In 
Lemma 6.1, we recall the basic equations of structure of [3] and some other prelim- 
inary results. In Lemma 6.2, we compute certain integral formulas modulo lower 
order terms. In Lemma 6.3, we give some variational formulas. Additional formu- 
las useful in studying manifolds with non totally geodesic boundaries are available 
from the authors upon request as are the proofs are entirely elementary if somewhat 
tedious. 

Lemma 6.1. 

1) dvolm = g 1 / 2 dx 1 ■ ■ ■ dx mi and Y vll a = \g ai {d^g vl + d u g m - d 1 g v , J ). 

2) IfD = -(g^d^ + P°d a + Q), then w v = \g v ^ + sTT^), 

E = Q- g v ^{d ll uj v + u„u>n - u cr T VIJ a ), fl ufi = d^u^ - d^u>, y + uj^lu^ - lo^lu^. 

3) L a[) = {d m ,V da d ) = T a[ r = d m g al3 . 

4) Pain —~ -^bb-.a ~~ ^ab-.bj Pam:rn Ramrnc;c H~ Pmrn-.ai Rabab\m — ^Rabrnb-.a^ Pij'.j — 
2^";*' Parn;a ~2^\m Pmm;m-! ^Pmm-.m — T\m Rabba\mi Rabcm — ^bc-.a ~~ ^ac\b-, 

5) f,a — f:ai f:am f;ma f-.m-.a^-^abf-.bi f:ab f:ab ^abf-.mi f:ii ~~ f:mm~^~f:aa 
^aaf-.m-) f-,aam — f\maa ~ Pamf-.a PmmJ ;m; f\maa — f\m:aa H~ ^ab:af :b H~ ^--^acf-.ac 
Lac^acf-jm ^aaf \mrrn f;iim — f;mmm H~ f;m:aa H~ ^--^ab-.af-.b ~~ ^aa-.bf-.b H~ ^--^acf-.ac ~~ 
Lac-Lac f;m Lddf-^^ pmm f\ m j f\aamm — f;rnrnaa^^Ramrnc:af:c ^Pmmj ^mm H~ 
^Rammbf-.ab ~ f;aPmm;a :brn Pmm;mJ;m- 



We must integrate by parts in computing certain integral invariants. The fol- 
lowing formulas are useful; we suppress lower order terms of length greater than 
3. 

Lemma 6.2. 

1) J f-.m-.aaF-.m — J Q^[{f;aam-P-,m F\m J -.a^bb-.a ^^mf-,aLab:b)- 

2) J j ";am^\am — I dM^ f\aam^\m ^-m f-a^bb-.a ~t~ ^-.m f-.a-^ab-.b •••)■ 
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4) f qm f;aamm — J Qj^{^f;m-^ab:ab f-.mmm-Laa 2^ m f :m "•)' 

We consider the variation 

D(e) = e- 2 *fD, g(e) = e 2e fg, F(e) = e~ 2 ^ F. 

To keep the boundary conditions constant, we set 

S(e) = e-^{Lu m (0)-Lu m (e) + S}. 

Let e, be an orthonormal frame for the tangent and cotangent bundles of M with 
respect to the reference metric g(0). Let ei(e) — e~ e) 'ei and e l (e) — e e) V be the 
corresponding frames for the metric 5(e). We remark that contraction and differen- 
tiation do not commute; 

£| £ =o(*ii) = (£| £ =o*)«-2/*« 
for example. Although T is not tensorial, its variation is tensorial. Let 

we keep the distinction between lower and upper indices since 
(4- e U=oT)ij k ^ (feU=o r hk; (we shall not need t ijk ). 



Lemma 6.3. 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 

19 
20 
21 

22 
23 
24 
25 
26 



{j^\e=0^)ij k — Sikf-j + Sjkf;i — Sijf- 



k- 



(£U= T) mm a = -/;„, and (£U= r) mm m = f. 

(■^\ e =oL)ab = —Sabf:m + fL a b- 



-o(E) 



-fS+\{m-2)f, m . 
-2/£+i(m-2)/ ;! 



del 

(j^\e=oR)ijke — Sikf-je + Sjef-ik — Suf-,jk ~ fijkf-,u + 2fRijki- 

0-R)mabm — 3abf;mm f,ab 4" 2 j ' Rmabm- 

{■^\ e =op)ij = (2 - m)f-ij - g tJ f,kk- 

e=o(Pmm) ^fPrnm f:aa H~ (1 TTl) f;mm- 

£ =o(r) = -2/T + 2(l-m)/ iH . 

e=o(-E;m) = —3fE;m + ~ 2)f;iim — 2/ ;TO -E. 

£ =o(T;m) = -3/r ;m 4- 2(1 - m)f. iim - 2/ ;m r. 
e=o(</> ; ii) = -2f<P ; u 4- (^| £ =O0); l4 4- (to - 2)f ;i (t) ;l . 

£=o(0;aa) = ~2f(j)-aa + (£|e=O0);aa + (to - l)/ ;m </> ;m 4" (m - 3)/ ;a ;a . 
£=o( ( / ) ;mm) = ~ 2/rf3 ;mm 4~ (^|e=O0);mm _ /;m0;m 4- / ;a 0;a- 

e— o(/-'mm;m) — Pmm-,m 4~ (1 TO) i /. mmm f-aam 4~ 2 f - a p arn 2f- m p rnrn . 

e— o(-^;mm) ^f^;mm ^f;mm^ ^ f;m^:m 4~" f-a^;a- 

e— o(-^;mmm) ^f-^;mmm 9 f-,m-^-,mm 4~ 3f; a F- am 7 f-mm-^;m 4~ f-^amE-a 

rnrnm-L • 

£ =o(-B;oo) = -4/S„„ + |(m-2)/ ;iioo -2/ ;oo f;+(m-7)/ i0 f; i0 + (m-l)/ irn £7 ;ro . 
£ =o(^:aa) = -4/£ :aa + ± (to - 2)/ ;J4:aa - 2/ :aQ £ + (m - 7)/ :a £ :a 

e— o(^;mm) : 4/_E' ;mrn + ^ 2) f^Hmm ^f\mm^ H~ /;a^;a 5f;m$;m- 

c=o(E, H ) = —AfE-u + \{m-2)f, ujj - 2f, ll E+{m-Q)f, a E, a + {m-6)f. m E. m . 

e=o(T-aa) = ~^fT- aa + 2(1 - m)f, iiaa - 2f- aa T + (l7l - 7)/ ;0 T ;a 4" (TO - l)/ ;m T ;m . 

£=o(^:qq) = -4/r :aa 4- 2(1 - m)f,ii. aa - 2f. aa r + (to - 7)/ :a r :a . 

e=o(^";mm) 4/T- mm 4~ 2(1 Wl)f:iimm 2/- mm T 4~ f;a^~:a 5/ ;m T- m . 

£ =o(T;ri) = -4/r ;ii 4- 2(1 - m)f;iijj - 2/ ;ii r 4- (to - 6)/ ;i r ;i . 



rfr 

iL 

rfr 
_d_ 
rfr 
_rf_ 
rfr 
_rf_ 
rfr 
_rf_ 
rfr 
_rf_ 
rfr 
_rf_ 
rfr 
_d_ 
rfe 

2/: 

_rf_ 

rfr 

_rf_ 

rfr 

_rf_ 

rfr 

_rf_ 

rfr 

_rf_ 

rfr 

_rf_ 

rfr 
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27) 

28) 

29) 

30) 
31) 
32) 
33) 
34) 
35) 
36) 
37) 
38) 

39) 
40) 
41) 
42) 
43) 
44) 
45) 
46) 
47) 
48) 

49) 
50) 
51) 

52) 
53) 
54) 
55) 

56) 

57) 

58) 

59) 
60) 



';mPam 
d I 



d , 

f;aPm,7 



3iU=o 
d I 

3iU=o 

d I 
d I 

d I 

diU=o 

diU=o 
d I 

diU=o 
d I 

d?U=o 
d I 

d?U=o 
(2 -to 
d I 

3iU=o 
d I 

3iU=o 
d I 

3iU=o 
d I 

d I 

d?U=o 
d I 

d I 

diU=o 
d I 

d^U=o 
d I 

diU=o 

d I 

diU=o 
5 (m- 

d I 

3iU=o 
d I 

3iU=o 

d I 

3iU=o 
i(m- 
d I 

d?l«=o 
d I 

diU=o 
d I 

diU=o 
d I 

diU=o 



) = ~4fp 

:a + (m- l)f- a P 



mm:aa 
'mm; i 



;aamm 

9mm;m 



(1 - m)/ ; 



mmmm 



2/; 



mrnHmm 



+ 



J mm • 
,J mm • 



+ (l-m)/ ; 

mmaa f]bbaa ^f;amPam ^f\aaPmm 
a ~\- ( jJl 1 ) ,/"; m Pmrn ; m • 
f J mm:aa mm:aa J ;aa: 66 + (l - m)/ 

:mm:aa ^ f :aaPmm~^~ 
(m - l)!:aPmm:a- 
Pmm:mm ) 

,;a ~t~ \aPam\m + 2/ 

tE) = -4/t£ + ±(to - 2)/,«t + 2(1 - 
£ 2 ) - -4/£ 2 + (m - 2)/ ;ii £. 
r 2 ) = -4./Y 2 + 4(1 - m)f. U T. 
p 2 ) = -4fp 2 - 2f. u r + 2(2 - m)f, ijPij . 

if, \jn- 8f., ijPij . 

PmmE) = -4/ 'p mm E - f. aa E + (1 - m)f- mm E + \ {m - 2)f-up v 
p mm r) = -4fp mm T - f. aa r + (1 - m)f- mm T + 2(1 - m)f. u p n 

RammbPab) — AfFLammbPab f;abPab $ \rnmT f:aaPmm~\~ 
) Rammb f;ab • 

Pmm) = Pmm ~ 2f;a,aPmm + 2(1 — Tn)f- mm p mm . 

RammbRammb ) = -4/J2 am in ammo - 2f, ab R 

ammb ^f-.mmpmm' 

^amQam) ' ^f^am^am- 

E. , m S) = -4fE, m S + ±(to - 2)f, iim S - 2f, m ES + \{m - 2)f. m E. m . 
ES 2 ) = -AfES 2 + ±(to - 2)f. ti S 2 + (to - 2)f. m ES. 
S 4 ) = -4/S 4 + 2(m - 2)/ ;m S 3 . 

t^S*) = -4fT- m S + 2(1 - m)f-u m S + ±(to - 2)/ ;m r ;m - 2f. m rS. 

PmmS 2 ) = --ifp mm S 2 - f-aaS 2 + (1 - m)f. mm S 2 + (m - 2)f. m p mm S. 

tS 2 ) = -4/tS 2 + 2(1 - TO)/ :ri S 2 + (to - 2)/ ;m rS. 
S : „aS) = -4/5 :aa 5 - / :aa 5 2 + (to - 5)/ ;a 5 :a 5 + |(m - 2)/ ;m:aa 5+ 

2)5':aa e /'; m* 

S :a S:a) = S-.aS-.a ~ ^f;aS :a S + (m — 2)/ ;m:a 5' :a . 

F\m-^;m) ^f^m^^m ^f^m^^^m 2 ^) f',iim^;m ^ f \m^\m^ • 

F- m ES) = -6fF. m ES - 2f, m FES + \{m - 2)f, u F, m S+ 

F. m S 3 ) = -6/F. m> S 3 - 2/. m F5 3 + § (to - 2)/. mJ F ;m 5 2 . 

^m^m) — 6 t /*-F; m T; m Qif-^FT-ui ~\~ 2(1 777^) f;iim^;m C ^ , f\m^\m^'- 

F. m rS) = -&fF, mT S -2f, m FrS + 2{l-m)f, li F, m S +\{m~2)f, m F. m T. 



F m Pm, m 5* ) — 6 / ' F- : m Pmm $ 2 m F p rn m S f;aaF^m$~^~ 



(1 - m)f- mm F- m S + 5(771 - 2)f. m F. m p n 

F-mS-.aa) — ~ 6/-F ;m S' :aa — 2f^ m FS M 

2)F- m f- m:aa . 



diU=o 
2 (to 

d I 

2 (to 

d I 

diU=o 

(to — 1 

d I 
d I 

f:aaF :rnm -\- (1 TTl) f;mm^;mi " 



S+{m- 5 ;/:„/•:„.. S : „ + 



F-mmE) — -6fF- mm E + ,f- a F- a E — 2f- mm FE — 5f- m F- m E+ 

2) f;iiF- nlrn . 

F-mmS 2 ) — —6fF [mm S 2 + f- a F- a S 2 — 2f. mm FS 2 — 5f- m F- m S 2 + 

) f;mF- : mmS . 

F- mm r) = -6fF. mm T+f. a F. a T-2f. mm FT-5f [m F. m T+2(l-m)f. ii F. r 



F. 



mmPmm 



) 6/ F. rnrn p mrn 

+ /; 

a" \aPmm run i F P in m -5/ \m^\mpT\ 
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61) ^ e |e— o(^;mmm^) ^fF^ mrnm S-\-3f- a F- am S 7 f-,mmF-mS 2f- mrnm FS 
Qf;mF-,mmS ~F f;amF- : aS -\- 2 [ffl 2) f-ynF-mmm' 

62) ^|e^o(^;mmmm) 6/-/^ mmmm Qf:mmmF- : m + f;mmaF- a ^-6f:mmF-,mrn + 
^f;amF- ; am 14/;m-^;mmm ~t~ 6/ ;a -F'- mma + 6/ ;a i? amm fr-F';fr 2f;mmmmF. 

63) ^| £ =o(i aQ T :m ) = -4fL aa T. m + (1 - m)f- m T- m + 2(1 - m)f. iim L aa - 2f. m L aa T. 

64) ^| £=0 (i aa:& ^) = -4/L aa:66 S + (l-m)/ ;ro:aa 5+(m-5)/ ;6 L aa:6 5-/ :66 L ao 5 + 
\{m~ 2)f :m L aa:bb . 

65) ^| £= o(ia6:cia6:c) = ~4/L a fc :c _L a b :c — 2f- m:c L aa:c + 4/ :e L e t,_L a b :a — 2/ ;c L a hL a h :c — 

66) ^le^oC^m^aaibb) = — QfF- m L aa: bb — 2f :m FL aa: bb + (1 — f- m -.bbF :m + 
(m - •">i/./,/ / : „, /.„„.,, - f:aaF- m L b b- 

67) |:|e = o(f;miat:nl)) = — 6/F ;m i a h :a fc — 2f. m FL a b :a b — f;m:bbF- m + (m — 2)f :e bF- m L e b 
~\~2{jYl 3) f:eF :m L ae:a f-.bbF-m^aa f:eF : mF a a:e ■ 

68) ^le-ot^mmm-^aa) 6/i^ mmm Z/ aa 9/- m -F' ;7nm Z/ aa + 3f-bF-b m F a a 
7 f-^mmF-mLaa ~t~ f;amF^ a Lbb 2f- mnlrn FL aa -\- (1 Wl) f ;mF;mmm ■ 

Remark. Branson and Gilkey [3, Theorem 7.2] also computed a„ for mixed bound- 
ary conditions. Their formula for 04 gave incorrectly the values /33 = —42 and 
/?4 = 6. Vassilevich [13] showed the correct values were = —12 and /3 4 = —24. 
The error in Branson-Gilkey arose from incorrectly applying the variational formula 
contained in Lemma 6.3 (4), which is valid for Neumann boundary conditions, to 
the more general context of mixed boundary conditions. If one were to consider 
mixed boundary conditions, the relevant variational formula for S would become 
3j| e= o(S') = —fS + (m — 2)/ ;m II + /2 and the other variational formulas involving 
S in Lemma 6.3 would then need to be changed accordingly. We wish to caution 
the reader concerning this point. 
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